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Abstract 

A method based on infinite parameter conservation laws is described to factor linear dif- 
ferential operators out of nonlinear partial differential equations (PDEs) or out of differential 
consequences of nonlinear PDEs. This includes a complete linearization to an equivalent lin- 
ear PDE (-system) if that is possible. Infinite parameter conservation laws can be computed, 
for example, with the computer algebra package ConLaw. 

MSC2000 numbers: 70.S10, 35.04, 35.A30 

1 Introduction 

With the availability of computer algebra programs for the automatic computation of all con- 
servation laws up to a given differential order of the integrating factors (as described in [Jj, |Hj) 
conservation laws have been found that involve arbitrary functions, i.e. infinitely many parame- 
ters. In this paper we show how based on such conservation laws a linear differential operator can 
be factored out of a combination of the original nonlinear partial differential equations (PDEs) 
and their differential consequences. Possible outcomes include 

• a complete linearization into an equivalent linear system, 

• a partial linearization in the sense that a linear differential operator is factored out, splitting 
the problem into a linear one plus a non-linear problem of lower order and often fewer 
independent variables (e.g. ordinary differential equations (ODEs)), 

• the derivation of at least one linear equation from a nonlinear system (with the possibility 
of deriving further linear equations for the new mixed linear-nonlinear system). 



I 



An advantage of the procedure to be presented is that conservation laws need not be given ex- 
plicitly in terms of the arbitrary functions. It is enough to have the conservation law determining 
conditions solved up to the solution of a system of consistent and necessarily linear PDEs which 
have arbitrary functions in their general solution. 

The content of the paper is as follows. After comments are made on the computation of 
conservation laws in section |2l the four computational steps of factoring out linear differential 
operators are illustrated using the Liouville equation in section^ Sufficient conditions for complete 
or partial linearizations are listed in section followed by a discussion of computational aspects in 
sectionlHl A generalization involving the introduction of potentials in terms of which a linearization 
becomes possible is explained in section In later sections IHl El an illustration is given of how 
the method works when nonlinear equations linearize to inhomogeneous equations or to triangular 
linear systems. Further examples where a complete or at least a partial linearization is possible 
are given in the appendix. 

In this contribution we concentrate on computational aspects of the method and give exam- 
ples for all of the above scenarios. An extension of the method discussing complete and partial 
linearizability through point and contact transformations will appear in a future publication PQ, 
with numerous new examples and a comparison with other linearization methods found in the 
literature. 

2 Notation 

We follow the notation in [3] and denote the original nonlinear partial differential equations as 
= A a , the dependent variables by u 13 , a, (3 = 1 . . . q and the independent variables by x\ i = 
I... p. In examples dealing with functions u = u(x,t) or u = u(x,y), partial derivatives are 
written as subscripts like u xy = d 2 u/(dxdy). If a formula already contains subscripts then <9, 
will be used for d/dx l . The multi indices j,k denote multiple partial derivatives like Uj which 
in our notation include u a . With j^J we denote the differential order, i.e. number of partial 
derivatives represented by j. Total derivatives with respect to x % will be denoted as Di. We apply 
the convention that summation is performed over terms that involve two identical indices, one 
subscript and one superscript. For example, the divergence of a vector field P % would be denoted 
as DiP 1 (= J2i DiP 1 ). The procedure to be presented repeatedly uses adjoint differential operators 
as follows. For given functions f A (x l ), A = l..r, let linear differential expressions Hk be defined 
as 

H k = a J kA djf A , k = l,...,s, 

with coefficients a{ A = (V) and summation over A and the multi index J. The corresponding 
adjoint operators H Ak are computed for arbitrary functions G k (x t ) by repeatedly reversing the 
product rule of differentiation for the sum G k Hk to get 

G k H k = ^HX^ + D.P 1 (1) 

where 

Hl k G k = (-l)* J dj (a J kA G k ) . (2) 
and P % are expressions resulting from integration by parts with respect to dj in this computation. 
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3 Conservation laws with arbitrary functions 



Conservation laws can be formulated in different ways (see jS] for four different approaches to 
compute conservation laws). The form to be used in this paper is 

AP ! = Q a A Q (3) 

where the components P l of the conserved current and the so-called characteristic functions 
Q a are differential expressions involving x l ,Uj. Other forms of conservation laws can easily be 
transformed into ©• One approach to find conservation laws for a given system of differential 
equations = A a is to specify a maximum differential order m of derivatives u 1 } on which P l , Q a 
may depend and then to solve condition (j3J) identically in x l , u 1 } for the unknown functions P l , Q a . 
Due to the chain rule of differentiation in Q the total derivatives Di introduce extra derivatives 
u% with j^K = m + 1 > m, i.e. derivatives not occuring as variables in P l ,Q a . Splitting with 
respect to these results in an overdetermined and linear system of PDEs for P l , Q a } 

What is important in the context of this paper is that a differential Grobner basis can be 
computed algorithmically and from it the dimension of the solution space can be determined, 
i.e. how many arbitrary functions of how many variables the general solution for P l , Q a depends 
on. In extending the capability of a program in solving conditions Q by not only computing a 
differential Grobner basis (for linear systems) but also integrating exact PDEs and splitting PDEs 
with respect to only explicitly occuring Uj (which here act as independent variables), the situation 
does not change qualitatively. The result is still either the explicit general solution or a linear 
system of unsolved PDEs 

Q = C k (a*,u^f A ), k = l,...,r, (4) 

for some functions f A (x^,Uj) where this system is a differential Grobner basis and allows one to 
determine algorithmically the size of the solution space. The functions f A are either the P\Q a 
themselves or are functions arising when integrating the conservation law conditions Q- 

If the conservation law condition Q is solved, i.e. P\ Q a are determined in terms of x\ Uj, f£ 
possibly up to the solution of remaining conditions (@J) then it is no problem to use a simple 
division algorithm to determine coefficients L k satisfying 

Q a A a = D t P* + L k C k (5) 

identically in x l ,Uj,f A . The coefficients L k are necessarily free of f A because Q is linear and 
homogeneous in Q a , P l and this property is preserved in solving these conditions, so C k are linear 
and homogeneous in ff as well and L k must therefore be free of ff . We will call relation © a 
'conservation law identity' because it is satisfied identically in all x\uj and ff. 

4 The Procedure 

The individual steps of our method are shown in detail to demonstrate that all steps are algo- 
rithmic and can be performed by computer. The Reduce package ConLaw has the algorithm 

^^Note that regarding © as an algebraic system for unknowns Q a implies division through A Q and does therefore 
not produce Q a which are regular for solutions u a of the original system A a = 0. For details regarding the ansatz 
for Q a see 0. 
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implemented and performs it whenever a conservation law computation results in a solution in- 
volving arbitrary functions possibly up to the solution of a linear system (J3J). 

Input 

Input to the procedure is the conservation law identity (J5J) 

Q a A a = D t P l + L k C k (6) 

including expressions for all its constituents Q a , P l , L k , C k in terms of x l ,u a , f A . 

To start the procedure the functions f A have to depend only on the variables x % . If they depend 
on Uj then a linearizatin will necessarily involve a change of variables. This case is treated in pQ. 

Step 1. 

If all functions f A depend exactly on all p independent variables x % then proceed to step 2. Step 
1 is concerned with the case that not all f A = f (x l ) depend on all x\ To add the dependence 
of, say f B on x J one has to 

• compute 

Z := (Q a A a - DiP* - L k C k ) 
which vanishes modulo = djf B and therefore must have the form 

Z = M J dj (f B ) 

with suitable coefficients M J and summation over the multi index J, 

• compute the adjoint Z* B as in (UJ),© to bring Z into the form 

Z = D i P i + d j f B Z* B , (7) 

• rename P % + P % — > P l and adds a new condition C r+ \ = djf B and multiplier L r+1 = Z* B 
to arrive at a new version of the conservation law identity Q a A a = DiP 1 + L k C k where the 
function f B depends now on xK 

This process is repeated until all f A depend on all x\ 

Example 1: We will illustrate the steps of the procedure with an investigation of the Liouville 
equation 

= A := u xy - e u . (8) 

Although it is not completely linearizable, we choose this equation because it involves computa- 
tions in each of the first three steps. 

For the Liouville equation a conservation law identity involving an arbitrary function f(x) is 
given through 

(f x + fu x )A = D x (-fe u )+D y (f x u x + ful/2), (9) 

i.e. Q = f x + fu x , P x = -fe u , P y = f x u x + fu 2 x /2, C k = 0. Adding a y-dependence to / 
requires to add to the right hand side of our identity (0) the terms 

Z fxyU x fy^x 
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which in adjoint form (J7J) read 

Z = D x (-fyU X ) + (U XX - U x 2 /2)fy 

giving the new conservation law identity 

(/* + fu x )A = D x (-fe u - f y u x ) + D y (f x u x + fu 2 x /2) + (u xx - u x 2 /2)f y . (10) 

Step 2. 

As the Q a are linear homogeneous differential expressions for the f A one can compute adjoint 
operators Q°£ as in (HJ),© by expressing 

Q a A a = f A Q a A *A a + D,P\ 

After renaming P % — P l — > P % the conservation law identity takes the new form 

f A Q a A *A a = D l P l + L k C k . (11) 

In the case of the Liouville equation we partially integrate 

(/, + fu x )A = f(u x - D X )A + D x (fA) 

and get the conservation law identity 

f(u x - D X )A = D x (-fe u -f y u x -fA) + D y (f x u x + fu 2 x /2) + (u xx -u x 2 /2)f y 

= D x (-f y u x - fu xy ) + D y (f x u x + fu 2 x /2) + (u xx - u 2 /2)f y . (12) 

Step 3. 

Because the are linear homogeneous differential expressions in the f A we can compute the 
adjoint form of L k Ck as in (JTJ) , (|2J) by expressing 

L k C k = f A CX k L k + D,P\ 

After renaming P % + P % — > P % the conservation law identity takes the new form 

f A Q a A *A a = + f A C2 k L k . (13) 

In our example partial integration gives 

(u xx - u 2 /2)f y = D y ((u xx - u 2 /2)f) - f(u xx - u 2 /2) y 

and substituted into ()12j) the new conservation law identity 

f{u x - D X )A = D x (-f y u x - fu xy ) + Dy (j ' X U x + f u^. 1 2 + f(u xx - u 2 /2)^ - f(u xx - u 2 /2) y 
= -f(u xx -u 2 /2) y (14) 

after simplification. 
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Step 4. 

This step does not involve any computation, it merely completes the constructive proof how 
linearizations are achieved. 

By bringing f A CX k L k to the left hand side of the conservation law identity (fT3j) we get 

f A (QfA a - CX k L k ) = D t P l (15) 

which still is an identity for arbitrary functions u a , f . Applying the Euler operator with respect 
to / (for its definition see e.g. to the left hand side of (fTB]) gives the coefficient of f A and 

on the right hand side gives zero as it is a divergence, 2 i.e. we get 

QX&a = CX k L k identically in u a for all A. (16) 

The vanishing of DiP 1 on the right hand side of (fTf)) was therefore not accidental. 
For the Liouville equation the identity (fTB*j) takes the form 

{u x -D x )A = -D y L = with (17) 

L = u xx -u x 2 /2. (18) 

Integrating at first (fTTjl to L = L(x) leaves the Riccati ODE 

u xx - u 2 x /2 = L{x) (19) 

for u x to be solved, for example, through a linearizing transformation u(x,y) = —2\og(v(x,y)). 
Output 

The result of the procedure are expressions QT,CX k and L k . The relation 

CX k L k = (20) 

is a necessary condition which can be solved by first regarding L k as dependent variables and then 
solving 

L k = L k {u a j) (21) 

for u a = u a (x l ). The system (J2"U|) . (}2"T]) is a sufficient condition for the original system A a = if 
Q°£ is an invertible algebraic operator and it is a complete linearization if (J21|) is purely algebraic 
in u a . 



5 Scope of the procedure 

The degree to which the original system A Q = can be linearized depends on properties of the 
conservation law identity that has been computed: the number of functions f A and the number of 

2 To prove this statement without Euler operator we could choose the f A to be zero outside some region R such 
that an integral over a volume with boundary outside R will vanish using Gauss law on the right hand side of 
identity (|15|l as P l are linear homogeneous in the f A . Because the f A are arbitrary inside R the coefficients of the 
f A on the left hand side must vanish identically. 
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variables each f A depends on, the differential order of derivatives of f A with respect to x l ,Uj in 
Ck and in Q a . Some properties, like the size of the solution space of remaining conditions (J1J) are 
essentially independent of the extent to which these conditions are solved. Other criteria, like the 
number of functions f A and the number of their arguments does depend on the extent to which 
conditions (HJ) were solved. The strength of the procedure to be presented is to be able to handle 
a wide range of situations. 

The following is a list of five scenarios, sorted from most special, fully algorithmic (and most 
beneficial) to most general, not strictly algorithmic (and less beneficial). We refer to the compu- 
tational steps described in section HI as 'the procedure'. 

• If the following criteria are met: 

1. the size of the solution space of = A a is equal to the size of the solution space of 

o = c k , 

2. the conditions = Ck involve q functions f A , (equal to the number of functions u a in 
A a ) and all f A depend on p variables (equal to the number of variables u a depend on), 

3. the functions Q a expressed in terms of f A involve f A only algebraically, i.e. no deriva- 
tives of f A , 

4. functions f A do not depend on jet variables Uj, i.e. f A = f A (x l ), 

then the procedure will algorithmically provide a linearization of the system A a = 0. 
Example 2: The Burgers equation in the form 

= Ai := u t - u xx - uu x (22) 

for a function u(x,t) can not be linearized but in the potential form 

= A 2 := v t - v xx - v x 2 /2 (23) 

for v(x,t) a conservation law identity involving a function f(x,t) is given through 

fe v ' 2 A = D t (f2e v ' 2 ) + D x (f x 2e v ' 2 - fe v ' 2 v x ) + 2e v ' 2 (-f t - f xx ) (24) 

and the related linearization is 

L = 2e v ' 2 
e v/2 A = L t -L xx = 0. 

A proof that every non-linear PDE (-system) that is linearizable through point or contact 
transformations can be linearized this way will be given in pp. 

• If criteria 1,2,3 are satisfied but not 4 then a linearization is possible but at the price of 
a change of variables, which will be a point or contact transformation if it is invertible or 
otherwise it will be a non- invertible transformation depending on derivatives of u a . Further- 
more, in all such cases the transformation can be derived explicitly from the conservation 
law identity as will be shown in pp. 
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• If criterion 3 is not satisfied then the partially or completely linearized equations may only 
be a necessary but not a sufficient condition for A a = 0. 

• If criterion 1 is satisfied but not 2 then 

— if functions f A of fewer than p variables occur then one can add extra variable depen- 
dencies through step 1 of the procedure, 

— if more than q functions f A occur in = Ck or functions f A of more than p variables 
occur then one has to integrate more of the conditions = Ck in order to be able to 
linearize the original system completely (a full treatment of this case will be given in 

• If criterion 1 is not satisfied but the solution space of Ck involves at least one arbitrary 
function of one argument then the method will result in a differential expression for Uj 
which vanishes modulo = A a and factorizes into a linear differential operator acting on 
a non-linear differential expression. Typically this leads to a PDE for u a which is lower in 
differential order than A a for one of the x l . In example 1 in section H] and examples 8,9,10 
in the appendix an equation in one less variable results, i.e. an ODE. 

The algorithmic beauty of the procedure is that the above wide range of situations are covered by 
one and the same 4-step algorithm. 

The case that a non-local linearization exists in which the L k depend on integrals of u a is 
not covered directly as the computer algebra package ConLaw does not compute non-local con- 
servation laws. On the other hand single conservation laws (without parametric functions) can 
be used to introduce potentials such that the original system re-formulated in these potentials is 
linearizable. This approach has been successful in all 6 linearizable evolutionary systems found in 
jS]. Examples given in this paper are the system (|2*5j) . (|2l)|) in the section the system (jSJ), 
in section El and the system (|fiU|) in the appendix. 

6 Computational Aspects 

Given a nonlinear PDE system = A Q , what are possible computational hurdles to be overcome in 
oder to find a linearization? The method described in section 0] is algorithmic and does not pose a 
problem. The formulation of conservation law conditions (jHJ) and their analysis through computing 
a differential Grobner basis = Ck is algorithmic as well and could only become difficult because 
of a growing size of equations. 

A first computational challenge lies in the fact that for linearizable systems = A a the conser- 
vation law conditions ([3J) have a general solution involving arbitrary functions. It is well known 
that systems of equations with a large solution space are much harder to solve than systems with 
only few solutions or no solutions. To incorporate many solutions, algebraic Grobner bases for 
algebraic systems have to be of high degree and differential Grobner bases for differential systems 
have to be of sufficiently high differential order. As a consequence, the length of expressions en- 
countered during the Grobner basis computation is more likely to explode and exceed available 
resources. 



8 



The second challenge is to integrate a Grobner basis = C*. sufficiently often to meet criterion 
2 in section^ Because the general solution of the conservation law conditions involves arbitrary 
functions, any integrations to be done can only be integrations of PDEs, not of ODEs. 

The package Crack that is used to compute the examples in this paper differs from similar 
other programs (as listed in in that it has a number of modules addressing the above prob- 
lems. For example, the growth of expressions is lowered by a module for reducing the length of 
equations by replacing them through a suitable linear combination of equations as described in 
P]. Integrations are handled by a module that integrates exact PDEs, that is able to introduce 
potentials to integrate certain generalizations of exact PDEs and that determines monomial in- 
tegrating factors to achieve integration (see (ID])- A relatively new module applies syzygies that 
result as a by-product of a differential Grobner basis computation. This module allows to perform 
integrations more efficiently and to avoid a temporary explosion of the number of functions of 
integration generated in the process (see PI]). A module that integrates underdetermined lin- 
ear ODEs with non-constant coefficients is often useful in the last stages of the computation. A 
description of the algorithm and its implementation is in preparation. 



7 An example requiring the introduction of a potential 

The following example demonstrates that a linearization of a non-linear equation or system may 
only be possible if it is reformulated in terms of potentials which in turn might be found by 
studying conservation laws. 
Example 3: The system 

= Ai := ut — u xx — 2vuu x — 2(a + u 2 )v x — v 2 u 3 — bu 3 — auv 2 — cu (25) 
= A 2 := v t + v xx + 2uvv x + 2(b + v 2 )u x + u 2 v 3 + av 3 + bvu 2 + cv (26) 

with u = u(x, t), v = v(x, t) and constants a, b, c results as one of 15 cases of a class of generalized 
non-linear Schrodinger equations j3] . This system itself does not have conservation laws involving 
arbitrary functions but it has the zeroth order conservation law 

vAi + uA 2 = D t (uv) + D x {v x u — u x v + bu 2 — av 2 ) (27) 

which motivates the introduction of a function w(x,t) through 

w x = uv, (28) 
— w t = v x u — u x v + bu 2 — av 2 . 

The remaining system to be solved for r := u/v and w simplifies if we substitute 

w = ^\ogz (29) 

with z = z(x,t). This substitution is not essential for the following but it reduces the size of 
the resulting system for r(x,t), z(x,t) and eases memory requirements in the computation of 
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conservation laws of the resulting system A 3 , A4: 

= A 3 := 2rr t zl + r 2 x z 2 x + 2ar 2 r x z 2 x - 2br x z 2 x + 2r 2 z x z xxx 

-r 2 z 2 xx + 2ar 3 z x z xx + 2brz x z xx + Acr 2 z 2 x 
= A 4 := r x z x + rzt - ar 2 z x + bz x . (30) 

The program ConLaw finds a conservation law with integrating factors 

Q 4 = r- 5 / 2 z x ' 3 / 2 (-2z x r(f x r-f x )-r x z x (fr + f) + z xx r(fr-f)) 

involving two functions f(x,t), f(x,t) that have to satisfy the conditions 

= d := -f t + f xx + cf- 2af x 
= C 2 := ft + L + cf-2bf x . 

The conservation law identity takes the form 

Q 3 A 3 + Q 4 A 4 = D t P* + D X P X + L X C\ + L 2 C 2 (31) 
with some conserved current (P t ,P x ) and coefficients L l ,L 2 of Ci,C 2 

L 1 = L 2 = A^Zfr. (32) 

Derivatives f x , f x in Q 4 can be eliminated by adding total x-derivatives 

D x (r- 5 / 2 z x -^ 2 2z x r(fr - /)A 4 ) 
to the left hand site of the identity (|31|) and to D X P X . The modified form of the identity (}3*T]) is 
= z- 3 / 2 r-^ 2 (2z x r(fr-f)D x A,-2r x z x (fr-f)A 4 + (fr + f)A 3 ) 
= A (Ajz~Tr(rf - /)) 

+D X {2z x -^ 2 r^l 2 {-2f x z x r 2 - 2f x z x r + r x z x fr - r x z x f 

+z xx fr 2 + z xx fr + Az x afr 2 + Az x bfr)^ 

+L 1 C 1 + L 2 C 2 . 

Partial integration of L X C\ + L 2 C 2 until /, / appear purely algebraically makes necessarily P l = 
P x = 0. Because /, / are free we obtain the identities 

= r- 3/2 z x r 3/2 (A 3 + 2rz x D x A 4 - 2r x z x A 4 ) = L] + L xx + cL l + 2bL 2 x (33) 
= r- 5 / 2 z x - 3/2 (A 3 -2rz x D x A 4 + 2r x z x A 4 ) = —L 2 + L 2 XX + cL 2 + 2aL l x (34) 

completing the linearization. For any solution L X ,L 2 of (JHB|), (JSH), equations ((HSj) provide r, z^. 
With Zt from ()30|1 we get z as a line integral, w from ()29|) and m, v from r and equation ()28j) . 

In the following section the effect of our method on PDEs is investigated that linearize to 
inhomogeneous equations. 
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8 Inhomogeneous linear DEs 



If the general solution of conservation law determining equations involves a number of free con- 
stants or free functions then individual conservation laws are obtained by setting all but one to 
zero. The remaining terms are homogeneous in the surviving constant or function. The question 
arises whether our conservation law based method is suitable to find linearizations that lead to 
linear but inhomogeneous equations. 

Example 4 : For the (ad hoc constructed) equation 

= A := 2uu t + 2uu xx + 2u 2 + 1 (35) 

the conservation law identity 

(f x + ft) A = D t (f x u 2 + f t u 2 + f) + 

D x {-f xx u 2 + 2f x uu x - f t , x u 2 + 2f t uu x + f) + 

^ (ft,x f xxx ftxx ftt) 

involves functions f(x,t),f(x,t) and establishes a conservation law provided /, / satisfy 

= ft,x ~ fxxx ~ ftxx + ftt- (36) 

Our method gives the linear system 

= D X A = L tx + L xxx (37) 
= D t A = L tt + L txx . (38) 
L = u 2 



The system (J3~7j) . (|HH|) represents the x and t derivatives of the linear equation 

= L t + L xx + 1 (39) 

which is equivalent to equation (|35|) and is an inhomogeneous linear PDE. Although our lineariza- 
tion method does not quite reach (|39|) . it nevertheless provides L = u 2 as the new unknown 
function which makes it easy to get to the equivalent linear equation (}3*9"j) through a change of 
dependent variables in (j33j) or through an integration of (jHZj), (J3~%j) . 

The way how homogeneous consequences can be derived from an inhomogeneous relation is 
to divide the inhomogeneous relation through the inhomogeneity, i.e. to make the inhomogeneity 
equal 1 and then to differentiate with respect to all independent variables and to obtain a set of 
linear homogeneous conditions in the same way as equations (|3*7j). (J3~%j) are consequences of (J3TI|) . 
The application in the following section leads to an inhomogeneous linear PDE with non-constant 
inhomogeneity. 
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9 An example of a triangular linear system 



A generalization of complete linearizability of the whole PDE system in one step is the successive 
linearization of one equation at a time. 

Example 5: Assume a triangular system of equations, like the (ad hoc constructed) system 

= Ai := u t (40) 
= A 2 := vv t — uvv xx - uv 2 (41) 

with one equation ()4U|) involving only one function, say u = u(x,t), and this equation being 
linear or being linearizable and a second nonlinear equation being linear or linearizable in another 
function v = v(x,t). How can the method in section |U be used to recognize that such a system 
can be solved by solving successively only linear equations? 

In determining all conservation laws for this system with unknown functions v, u and with 
integrating factors of order zero we get apart from two individual conservation laws with pairs of 
integrating factors (Q l , Q 2 ) = (\, — -) and (^-, — — ) only one with a free function f(u, x): 

/«Ai = DJ 

which indicates the linearity of A x but not the linearity of A 2 in v once u(x,t) is known. 

The proper way of applying the method of section 0] is to compute conservation laws of = A 2 
alone which now is regarded as an equation for v(x, t) only. The function u(x, t) is assumed to be 
parametric and given. We obtain the identity 

2/A 2 = D x (f x uv 2 + fu x v 2 - 2fuvv x ) + D t (fv 2 ) - v 2 {f t + uf xx + 2u x f x + u xx f). 

which is a conservation law if / satisfies the linear condition 

= ft + ufxx + 2u x f x + u xx f. (42) 

This provides the linearization 

= 2A 2 = L t - uL xx 
L = v 2 . 

The reason that now a linearization of A 2 is reached is that in the second try u is assumed to be 
known and therefore u,u xx , . . . are not jet-variables and hence the condition ()42j) has solutions, 
otherwise not. 

Examples where this triangular linearization method is successful are the systems (17), (18) 
in [Sj. We demonstrate the method with one of these (system (17)), the other is similar. 
Example 6: The system 

= Ai := ut — u xx — Auvu x — Au 2 v x — 3vv x — 2u 3 v 2 — uv 3 — au (43) 
= A 2 := v t + v xx + 2v 2 u x + 2uv v x + 2u V + v 4 + av (44) 
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involves functions u(x,t),v(x,t) and the constant a. The single conservation law 

= vAi + uA 2 = D t {uv) + D x {uv x — u x v — u 2 v 2 — v 3 ) 

motivates the introduction of a function w(x,t) through 

w x = uv (45) 
— w t = uv x — u x v — u 2 v 2 — v 3 . (46) 

Substitution of u from (|43j) brings equations (jH)) and in the form 

= A 3 := w t - -(-2v x w x + w xx v + w 2 v + v 4 ) (47) 

v 

= A 4 := v t + v xx + 2w xx v + 2w 2 v + av + v 4 . (48) 

This system obeys conservation laws that involve a function f(x, t) that has to satisfy f t = f X x+af- 
Our procedure provides the linearization 

e w {vA 3 + A A ) = L\ + L l xx + aL l = (49) 
L 1 := ve w . 

The second linearized equation can be obtained by 

• substituting v = L 1 /e w into equations (}4Tj) and flHj) : to get the remaining condition 

= A 5 := w t - w xx - 3w 2 + 2w x L l x {L l Y l - (L 1 ) 3 e~ 3w , (50) 

• assuming L 1 has been solved from (f47Ij) and treating L 1 (x,t) as a parametric function when 
computing conservation laws for equation ()50|) which turn out to involve two functions that 
have to satisfy linear PDEs, 

• performing the linearization method to find that the remaining equation (|5T)|) linearizes with 

L 2 = e 3w to 

e 3 -A 5 = L 2 - L 2 XX + 2LILI/L 1 - 3(L 1 ) 3 . (51) 

Because the condition (JoTj) is inhomogeneous for L 2 due to the term 3(L 1 ) 3 , actually two homo- 
geneous linear equations are generated which are the x— and t— derivative of (joTj) divided by 
3(L X ) 3 (see previous section about linear inhomogeneous equations). But as the function L 2 = e 3w 
results in this process, it is no problem to find (J51j) from (J5U|) directly or from an integration of 
these two equations. This completes the linear triangularisation of the original problem (|4*3*j) . (j44j) 
to the new system (j4Uj) . (J5T]) . 
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10 Summary 



The paper starts with introducing conservation law identities as a natural way to formulate infinite 
parameter conservation laws. 

Conservation law identities are the input to a four step procedure that returns a differential 
consequence of the original system together with a linear differential operator that can be factored 
out. 

Sufficient conditions on the conservation law identity which either guarantee a complete lin- 
earization or at least a partial linearization are discussed. 

The possibility to find a non-local linearization arises from the application of single (finite 
parameter) conservation laws with the aim to introduce potentials which satisfy infinite parameter 
conservation laws and thus allow a linearization. 

In examples it is demonstrated how the standard procedure can lead to inhomogeneous linear 
PDEs and how a successive linearization of one equation at a time may be possible when the 
whole system can not be linearized at once. 



Appendix 

In the appendix we list further examples of linearizations and integrations without giving details 
of the calculations. 

The first example of the Kadomtsev-Petviasvili equation demonstrates what our method gives 
when a PDE has p independent variables and the conservation law involves free functions of less 
than p — 1 variables. Although the result will be less useful than in the other examples, we still 
include it for illustration. 

Example 7: The Kadomtsev-Petviasvili equation 

A Ui x -\- U X xxx ^U XX U -\- 2lt x . ^"yy 

for u(t,x,y) has four conservation laws with a zeroth order integrating factor and an arbitrary 
function f(t) as given in |B]. We comment on one of these four with an integrating factor fty 3 +Qfxy 
as the situation for the others is similar. Omitting the details we only give the result of our method: 

L 1 = y [utxxxV 2 + Zutxuy 2 + u a y 2 + 2u t u x y 2 

—QutX — Qu xxx x + Qu xx — Ylu x ux + Qu 2 ^j 

L 2 = —u ty y 3 + 3u t y 2 + 6u y xy — 6ux 
y(6xA - y 2 D t A) = L l x + L 2 y . 

The arbitrary function f(t) involves only one independent variable t and the conservation law 
= L\ + Ly involves two functions L l ,L 2 and has derivatives with respect to two variables x, y 
and is therefore not as useful as if it would be a single total derivative. 

The three following equations were shown to the author first by V. Sokolov [4j who obtained 
their integrations earlier and independently. We add them here to demonstrate that these results 
can be obtained in a straight forward procedure. 
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Example 8: For the equation 

= A:=u xy - e u ^/u x 2 - 4, u = u{x,y) (52) 

a conservation law with an arbitrary function f(x) enables to factor out D y leaving an ODE to 
solve 

(«. - D x ) ( = D y ( + = o. (53) 

1 lu 2 - 4 J \ Ju 2 - 4 J 



Another conservation law with an arbitrary function g(y) gives 

U r e U „ | , „ / 1 n 1 



Uy ~ /z^ ~ Dy } A = Dx {~~ Uyy + ^ u y + 2 e2u) l a 



Example 9: For the equation 

= A := - [ — 1 — ] u x u y , u = u(x, y) (55) 

\u — x u — y J 

a conservation law with an arbitrary function f(x) similarly to the above example provides 

y — x A „ /A\ ( u rr 2{u r — 1) u \ , . 

J -A + A,. I — ) = D„ — - - - — = 0. (56) 



(u — x)(u — y) \u x J \u x u — x (u — x)x / 

A second conservation law is obtained from an arbitrary function g(y) and is equivalent to (|56|) 
after swapping x «-> y. 

Example 10: For the equation 

2 

= A := u xy - - ^ y/u x u y , u = u(x, y) (57) 

a conservation law with an arbitrary function f(x) gives 

1 ; 1 1 \ . _ / A \ „ / u xx 2t/u x ~\ 



A + D x I — — \ = D y I —p=z H — = 0. (5 



A second conservation law is obtained from an arbitrary function g(y) and is equivalent to 
after swapping x «-> y. 

The final example shows a linearization of a system that resulted in classifying non-linear 
Schrodinger type systems in [5]. 

Example 11: The system 

= Ai := Ut — u xx — 2vu x — 2uv x — 2uv 2 — u 2 — au — bv — c (59) 
= A 2 := v t + v xx + 2vv x + u x (60) 
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involves functions u(x,t),v(x,t) and the constants a,b,c. The trivial conservation law 

= A 2 = D t (v) + D x (y x + u + v 2 ) 

motivates the introduction of a function w(x,t) through 

w x = v (61) 
— w t = v x + u + v 2 . (62) 

Substitution of u, v from (JBTj) and (|52*jl brings equation (JoUj) in the form 

= A 3 = w tt + w 2 - w t a - w xxxx - 4w xxx w x - 3w 2 xx - Qw xx w 2 x - w xx a 
—w x — w 2 a + w x b + c. 

This equation admits a conservation law identity 

fe w A 3 = D t i(e w )tf-e w f t -e w fa] + 

D x [ (c ) xxx f ~\~ (c )xxfx ifi )xfxx & f xxx 

~(ae w )xf + ae w f x + be w f] 
+e w [ftt + aft - f xxxx - af xx - bf x + cf] . 

From this follows the linearization 

e w A 3 = Ltt — Ltd — L xxxx — L xx a + L x b + Lc = 
L = e w . 
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